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ABSTRACT 

Although the relationship between the far-infrared and cm-wave radio luminosities of normal 
galaxies is one of the most striking correlations in astronomy, a solid understanding of its physi- 
cal basis is lacking. In one interpretation, the "calorimeter model," rapid synchrotron cooling of 
cosmic ray electrons is essential in reproducing the observed linear relationship. Observed radio 
spectra, however, are shallower than what is expected of cooled synchrotron emission. In Thomp- 
son et al. (2006), a simple parameterized model is presented to explain how relatively shallow 
observed spectra might arise even in the presence of rapid synchrotron cooling by accounting 
for ionization losses and other cooling mechanisms. During the commissioning of the 42-element 
Allen Telescope Array, we observed the starburst galaxies M82, NGC 253, and Arp 220 at fre- 
quencies ranging from 1 to 7 GHz, obtaining unprecedented broadband continuous radio spectra 
of these sources. We combine our observations with high-frequency data from the literature to 
separate the spectra into thermal and nonthermal components. The nonthermal components all 
steepen in the cm-wave regime and cannot be well-modeled as simple power laws. The model 
of Thompson et al. is consistent with our M82 results when plausible parameters are chosen, 
and our results in fact significantly shrink the space of allowed model parameters. The model 
is only marginally consistent with our NGC 253 data. Assuming the Thompson et al. model, a 
steep electron energy injection index of p = —2.5 is ruled out in M82 and NGC 253 to >99% 
confidence. We describe in detail the observing procedures, calibration methods, analysis, and 
consistency checks used for broadband spectral nbse rv3tinns with tlj^ Mie n Telescope Array. 
Introduction 



The correlation between the far-infrared (FIR) 
and centimeter-wave radio emission of normal 
galaxies has been reco gnized for almost four 
decades ( van der KruitI [l971). It holds over 
more than four orders of magnitude in LpiR 
(jYun et aZ.ll200lh . is constant out to at least z — 1 



(jAppleton et aLl 12004) . and is obser ved in a large 
variety o f galaxies ( CondonI 19921 ) . The near- 
linearity ( Yun et al. 200 ih and universality of the 
FIR-radio correlation strongly suggest that it is 
intimately linked to the fundamental astrophysics 
of galaxies. 

The physical underpinning of the correlation 
is broadly understood to be the formation of 
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massive stars, which can believably be associated 
with the dominant contributions to both the FIR 
and contin uum cm-wave radio emissi on of galax- 
ies (Harwit fc Pacinil 119751 : IWunderhch fc Klein 



19881) . The FIR contribution is due to the hot 



dust envelopes surrounding these stars, which are 
energized by the copious amounts of Lyman con- 
tinuum radiation that massive stars emit. The 
radio contribution includes two components. The 
first is synchrotron emission from relativistic cos- 
mic ray (CR) electrons spiraling in the interstellar 
magnetic field. Such electrons can be accelerated 
in the supernova (SN) shocks that occur when 
massive stars die. The second, less significant, 
component is free-free emission, which predomi- 
nantly comes from H II regions. These too are 
associated with massive stars. 

This simple picture, however, neglects many 
important details. The timescales and dynamics 
of the radio and FIR emission are very different: 
CR electrons radiate over a timescale of 10''-10^ yr 
(see below) , during which they diffuse through the 
galaxy, are possib ly advected out of it via a galac- 
tic wind ( Breit schwerdt et aZ.lll99lh . and cool via 
a variety of mechanisms, including inverse Comp- 
ton (IC) scattering, ionization, free-free emission, 
and adiabatic expansion. Meanwhile, the FIR 
emission is prompt and comprised of compara- 
ble warm (~40 K) and cool ("cirrus," ~20 K) 
dust components, the latter of which is primar- 
ily heated by the i nterstellar radiation field, not 
newly-formed stars ('Walterbos & Schwerin 3 119871: 
|j3oulangcr & Pcrault 1988; Fitt et al. 1988). 

Much theoretical and observational effort has 
been expended in the cause of und erstanding the 
FIR-r ad io correlation i n deta il, 
(Il985l) 



Helou et al. 
Hummel et_a l! ('1988') 
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(ll985^: ICox et all 
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Price fc DuriG (1992): Lisenfeld et al (1996)). A n 
influential model has been that of Volkl ( 19891 ) . 



which presents a theory for understanding the 
FIR-radio correlation in the case that CR electrons 
are accelerated only by SNe and cool completely 
within a galaxy. In this situation, the galaxy 
acts as a "calorimeter" with regard to the energy 
injected into CR electrons. For the FIR-radio cor- 
relation to hold, the synchrotron luminosity per 
unit energy injected into CR electrons must be a 
universal constant. Taking the primary cooling 
mechanisms to be IC scattering and synchrotron. 



it is found that Us/iUph + Ub) must be nearly 
a constant value, where Ub = (^^/Stt) is the 
mean energy density of the interstellar magnetic 
field and Uph is the mean energy density of the 
interstellar radiation field. 

Bec ause the synchrotron co oling time scales as 
5^3/2 ([Thompson et a/.ll2006l hereafter T06), the 
applicability of the calorimeter model depends on 
the strength of the interstellar magnetic field. In 
a stronger field, the calorimeter model is more 
likely to hold as CR electrons will cool more 
rapidly and have less chance to escape the galac- 
tic disk before radiating all of their endowed en- 
erg;^ This is precisely the scenario suggested in 



T+Od where it is argued that starforming galax- 
ies have much stronger magnetic fields than com- 
monly assumed. The work assesses the applica- 
bility of the calorimeter model by considering the 
timescales for synchrotron cooling and escape from 
the galactic disk by CR electrons. Using the usual 
assumptions, the two are c ompara ble: using the 
stronger magnetic field that T -1-061 argues for, the 
synchrotron cooling timescale is much shorter, fa- 
voring the assumptions of the calorimeter model. 
T-H06I specifically argues that in starburst 



galaxies the magnetic field is dynamically relevant: 
its energy density is comparable to the hydrostatic 
pressure of the ISM, or 
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where is the surface density of the galactic gas 
disk. This is in co ntrast to the classic "minimum - 
energy" estimate (Burbidge 19561 : Longair Il994l) . 
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where rj is the ratio of total CR energy to energy 
in the CR electrons alone, L^, is the synchrotron 
luminosity, V is the estimated volume of the emit- 
ting region, and v is the observing frequency. For 
the parameter values commonl y assoc iated with 
starburst galaxies, B^q ^ Smin ( T-l-061 . fig. 1). 

An impor tant ob servational constraint on the 
argument of IT-H06I is the shape of the galactic 
nonthermal spectrum. It is typically argued that 
galactic spectra are insufficiently steep to be com- 
patible with the calorimeter mo del: obse r vation s 
indicate a « -0.8 {S^ oc z/"; ICondonI (|l992l ): 



Niklas et al. i 1997t )) rather than a ~ —1 as would 
be expe cted from the com bination of SN shock 
physics (jBlandford fc Eich lcr 1987) and spectral 
steep en ing ca used by strong cooling ([Condon 



19921) 



contends that more detailed anal- 
ysis mitigates this issue. As CR electron energy 
decreases, the timescales of ionization losses and 
free-free emission are, respectively, shorter and 
unchanged, so both of these proces ses will tend to 
flatten observed spectra. In simple steady- 

state simulations were performed to quantify this 
effect (cf. its fig. 3 and §4.3). These simulations 
lead to the general prediction that the nonthermal 
radio spectra of starburst galaxies should steepen 
with increasing frequency, asymptoting to a < —1. 

Most spectral index measurements are based on 
only a handful of fluxes. To accurately measure 
the steepening of a nonthermal spectrum, many 
more data are needed, because both the steepening 
and the ther mal contribut ion must be deter mined. 
The work of iPuric et al. and Niklas et al. 



(|l997l ) shows that determination of the thermal 
fraction is difficult and requires high-quality high- 
frequency data, even if a simple power-law non - 
thermal spectrum is assumed. As ICond"onl ( 19921 ) 
points out, the thermal fraction and nonthermal 
steepening can also be highly degenerate (cf. his 
Figure 6). One must sample a large number of 
points over a wide frequency range in uniform con- 
ditions to eflectively constrain these parameters. 

Though it is still bei ng commissio n ed, th e Allen 
Telescope Array f ATA: IWelch et all (|2009f) ). with 
its continuous frequency coverage of 0.5 GHz to 
10 GHz, is well-suited to testing the predictions 
of lT-]-06l by observing high-resolution, continuous 
spectra of starburst galaxies. Section [2] describes 
observations of M82, NGC 253, and Arp 220 made 
with the ATA in the fall and winter of 2008. In g3] 
our data reduction methods are described along 
with results for standard calibrator sources. Sec- 
tion[5]describes our modeling of the observed spec- 
tra to extract thermal and nonthermal contribu- 
tions. Finally, in fj5]we discuss our results and in 
^we present our conclusions. 

2. Observations 

We observed the starforming galaxies M82, 
NGC 253, and Arp 220 with the Allen Telescope 
Array on seven nights over the course of 2008 



September to 2009 June. The ATA is a large-A^", 
small-D (LNSD) cm-wave array located in Hat 
Creek, California, and is a joint project of the 
SETU Institute of Mountain View, CA and the 
Radio Astronomy Laboratory of the University of 
California, Berkeley. The scientific motivators and 
technical capabiliti es of the ATA are described in 
Welch aii (|2009[ ). At the time of the obser- 
vations the ATA was undergoing commissioning, 
with many aspects of the array being actively de- 
veloped. Forty-two antennas and one correlator 
had been installed at the site. Two-thirds of the 
antennas were typically available for use during 
our observations. 

We used the continuous frequency coverage of 
the ATA to observe the galaxies at almost a decade 
of sky frequency in a sequence of snapshot observa- 
tions. The observing frequencies ranged between 
1.0 GHz and 7.0 GHz with a spacing of 100 MHz, 
skipping some frequencies due to known radiofre- 
quency interference (RFI) or low priority. The 
primary beam sizes of the ATA at 1.0 GHz and 
7.0 GHz are ~210' and ~30', respectively, and 
the synthesized beam sizes are typically ^250" 
and ~35". The snapshots lasted between 60 and 
130 seconds, with the correlator dump time always 
being 10 seconds. All observations used a band- 
width of 104.9 MHz and 1024 correlator channels 
and sampled both vertical and horizontal linear 
polarizations. Observations of the science targets 
were interleaved with observations of the calibra- 
tors 3C 48, 3C 147, and 3C 286. These calibra- 
tors are compact to the ATA and have up-to-date 
broadband flux density models available from the 
vlaI Because of the need to set the flux density 
scale at numerous frequencies and the desire to ob- 
tain a large amount of data for sanity-checking, a 
large portion of the observing time was dedicated 
to calibrator observations. The total integration 
time on the calibrators was ten hours, compared 
to nine hours spent on the target galaxies. An 
observing strategy that gathered fewer redundant 
data could spend approximately half as much time 
observing calibrators as observing science targets. 



^Search for Extraterrestrial Intelligence 

^http : //purl . oclc ■ org/net/pkgwpub/vla- analytic- cal-models | 
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2.1. Details of observing methods 

Because the observations presented herein use 
a novel technique on a novel instrument, we shall 
describe our methods in detail. Each night was di- 
vided into several segments during which a subset 
of the targets and calibrators would be observed 
as a group. The typical duration of a segment was 
a few hours. Over the course of each segment, 
an observing script would repeatedly slew to each 
of the selected sources in turn, observing each at 
several frequencies before moving on to the next, 
until the segment's allotted time had elapsed. This 
scheme attempted to balance the desire to obtain 
calibrator observations temporally close to the sci- 
ence observations with the desire to not spend too 
much time moving from one target to the next. 
Most segments included observations at 1.4 GHz 
and 5.0 GHz to enable checks of repeatability and 
stability. The observations presented in this paper 
are summarized in Table [T] 

Observing over a large range of frequencies re- 
quires not just a large number of calibrator scans 
but also reconfiguration of the ATA hardware. In 
particular, the sensitivity of the log-periodic ATA 
feeds is a function of position along the axis of 
their pyramidal structure, with the point of peak 
sensitivity approaching the apex of the pyramid 
with increasing frequency. To allow for maximal 
performance, the ATA feeds can be focused by 
moving them in a piston motion al ong the nor- 
mal t o each antenna's focal plane (see Welch et all 
20091 fig. 4). The optimal focus positions have 
been determined theoretically and calibrated so 
that an ATA observer can specify the desired focus 
position in frequency space. The focusing opera- 
tion, however, is slow. In our work, before each ob- 
servation we set the feed focus frequency to the ob- 
serving frequency rounded up to an integral num- 
ber of GHz {e.g., 4.0, 4.4, and 4.9 GHz all map to 5 
GHz). This scheme was chosen because the sensi- 
tivity of the feed to frequencies above the nominal 
focus frequency drops off more rapidly than the 
sensitivity to frequencies below it. Within each 
segment we chose observing frequencies that lay 
within a single GHz span to avoid superfluous re- 
focus operations. The focusing and other time- 
consuming operations are parallelized as much as 
possible to conserve observing time. 

Because neither the feed sensitivity nor the 



noise power entering the ATA are spectrally flat, 
when the focus position of a feed changes, the to- 
tal power entering the analog system also changes. 
Just before digitization, the downmixed and fll- 
tered analog signal is passed through an "attem- 
plifier" (attenuator/amplifier) stage to adj ust the 
incoming power for optimal quantization ([Backer 
20071 ) . The amount of attenuation or amplification 



must be explicitly configured, but the best settings 
can be calibrated automatically during observa- 
tions. Clearly, any two observations that are to 
be compared or used for fiux calibration must be 
made with the same attemplifier settings. Fur- 
thermore, the attemplifier calibration step is also 
somewhat slow (~60 s). In these observations, 
snapshots made at similar frequencies shared at- 
templifier settings. Specifically, one set of attem- 
plifier settings was used for each integral GHz span 
{e.g., 4 < I'GHz < 5). The first time any frequency 
in a span is observed, the attemplifier settings are 
calibrated for that frequency and saved; for subse- 
quent observations during the same segment, the 
saved settings are reapplied. 

As an example of all of these considerations, 
consider segment #12, as detailed in Table[5] This 
segment was scheduled for observation from 03:42 
UT to 05:12 on 2008 Oct 15. (All dates and 
times are in UT.) The selected sources were 3C 48, 
3C 147, and M82. The timing and (short) dura- 
tion of this segment were set by the visibility of 
the science targets to the ATA: Arp 220 set just 
before this segment began, and NGC 253 rose just 
after it ended. Because M82 is circumpolar at the 
ATA's latitude, the other two targets were prior- 
itized whenever they were visible. The selected 
observing frequencies were 1.4, 2.4, 2.5, 4.3, 4.7, 
and 5.0 GHz. As described above, the observing 
script slewed to the sources in turn and observed 
each at several frequencies. The frequencies were 
chosen and grouped so as to reduce the amount of 
time spent adjusting the feed focus positions and 
attemplifier settings. At 04:08, the script would 
have observed 3C 147 at the trio of lower frequen- 
cies if possible, but that source had not yet risen. 
Instead, the script skipped those observations and 
looped back to observe 3C 48 at the trio of higher 
frequencies. 

Including time lost due to slewing, hardware 
initialization, software overhead, and so on, we 
made ~18 snapshots per hour while observing, in- 
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tegrating for 38% of the 53 hours spent using tele- 
scope resources during the observations presented 
in this paper. The efficiency of future observations 
will be higher: we have made an unusually large 
number of calibrator observations and have made 
no effort to reduce slew time between sources. Fur- 
thermore, our observations were made before the 
installation of a second independently-tunable cor- 
relator at Hat Creek, which doubles the observa- 
tional throughput. 

3. Data Reduction 

The data were reduced with a combination o f 
the analysis package MIRIAeH <|Sault et oi.lli995h 
and custom Pythoifl scripts based on the NumPy 
and IPython packages. Over the course of this 
project and other work on the ATA, significant 
improvements have been made to MIRIAD to ease 
its compilation, improve its ability to handle large 
datasets, and fix various bugs. A special module 
was written to allow easy management of and ac- 
cess to MIRIAD datasets in the Python language^ 

The two orthogonally-polarized signals that are 
received by the ATA feed elements are amplified 
by separate low-noise amplifiers (LNAs) and sub- 
sequently treated as entirely separate signals. The 
ATA correlator at the time of the observations 
accepted 64 such antenna-polarization signals, or 
"antpols" for short, making it capable of corre- 
lating 32 antennas with full-polarization coverage, 
64 antennas with single-polarization coverage, or, 
as encountered in practice, an intermediate case. 
Each ATA antenna is identified by a small inte- 
ger and a single letter, e.g. 3g. (These specify the 
node and element identifier of the antenna within 
the ATA architecture.) An antpol is specified by 
suffixing the antenna name with the letter X or 
Y, e.g. 2kY. We refer to a pair of antpols as a 
"basepol." 

Although the ATA correlator produces cross- 
hand polarization data, the calibration and anal- 
ysis of these data have not yet been investigated 
and tested at the time of this work. In this paper, 
cross-hand data were discarded, and all sources 
were assumed to be completely unpolarized. All 
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images produced were of the Stokes I parame- 
ter. Because the linear polarization of 3C 286 is 
~10%, this treatment will lead to ~1% errors in 
some measurements. This effect contributes to, 
but does not dominate, the typical uncertainty in 
our measurements (see ij3.4p . 

3.1. Flagging 

The ATA does not yet have an online flagging 
system and our observations were performed out- 
side of the protected bands. Fairly extensive man- 
ual editing of the data was performed. Several 
automatic criteria were applied. The 100 channels 
on both ends of the 1024-channel spectral window 
were flagged because the ATA digital filter aliases 
in large amounts of noise at the bandpass edges. 
Shadowed antennas were flagged, taking the effec- 
tive antenna radius to be 7.5 m (somewhat more 
than the 6.1 m physical diameter of the primaries). 
Finally, visibilities with ^/u^ + v"^ < lOOA were 
flagged as these short baselines were consistently 
sensitive to terrestrial and solar RFI, undesired 
large-scale structure, and possible cross-talk. 

A fraction of the antpols had problems with 
their cryogenics, LNAs, or delay calibrations, re- 
sulting in data for the night that could easily be 
rejected after a cursory look. Occasional correla- 
tor hardware failures resulted in wildly erroneous 
amplitudes in a small subset of the spectral data, 
e.g. the first quarter of the spectrum of all visibil- 
ities involving antpols laX, IcY, and 4kX. Such 
data were also flagged and the issue reported. 

After an initial check for bad antpols or cor- 
relator hardware, RFI was flagged with a semi- 
automated, channel-based routine. All of the am- 
plitude spectra at a given frequency observed over 
the course of a single segment were summed, re- 
gardless of time, baseline, polarization, or point- 
ing. A simple peak-finding routine identified sus- 
pected RFI channels by searching for extremely 
high amplitudes. The suspect channels were pre- 
sented to the user graphically and made available 
for interactive editing before a list of channels to 
flag was recorded. The listed channels for each 
frequency were flagged out for all observations at 
that frequency in the segment. Figure [T] shows an 
example of the summing and peaking-finding pro- 
cess. Figure [5] characterizes the RFI environment 
of Hat Creek via a schematic representation of the 
summed amplitudes across the entire observed fre- 
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quency range. 

This procedure is obviously suboptimal in the 
case of time-varying RFI, but it was found that 
most of the interference was virtually always 
present. Informal tests of more sensitive method- 
ologies that attempted to better mitigate time- 
variable RFI failed to improve results and were 
much more time-intensive to the user. 

A few bands had nearly omnipresent RFI that 
made astronomical observations nearly impossi- 
ble. These were: 500-1000 MHz (terrestrial televi- 
sion, radar, pagers, aeronautical navigation, self- 
interference), portions of L band (1150-1280 MHz, 
aeronautical navigation and GPS L5; 1520-1610 
MHz, which includes GPS LI at 1575 MHz), the 
lower end of S band (satellite communications, 
2150-2350 MHz, especially Sirius XM radio cen- 
tered at 2339 MHz), and C band (satellite down- 
links, 3700-4200 MHz). 

After RFI-afflicted channels were flagged, vari- 
ous manual data-quality checks were run to check 
for remaining issues. These included phase-closure 
rankings of the calibrator observations, evaluation 
of the amplitude and phase variability of the base- 
lines, and by-eye assessments of spectra. Many 
of these checks resulted in a few antpols or base- 
pols being flagged for the duration of the segment, 
while occasionally a temporary system glitch was 
revealed and edited out. 

Finally, after system temperatures were com- 
puted (see below), antpols with a Tgys above 300 
K were flagged. At the end of the flagging pro- 
cess, ^--^40% of the antpols and ~75% of the visi- 
bility channels in a typical dataset were discarded. 
Figure [3] shows the average data retention rates as 
a function of frequency. The deleterious effect of 
RFI is apparent, as is the moderate degradation 
of system performance with increasing frequency. 
Approximately 30% fewer antpols and visibilities 
are retained at 7 GHz as compared to 1.4 GHz. 

Ongoing retrofits to the array are increasing 
the amount of usable data generated during each 
observation. A large portion of this work com- 
prises improved RFI mitigation techniques, in- 
cluding better shielding, especially of the signal 
processing room; improved techniques and algo- 
rithms for interference excision; and stricter on- 
site RFI protocols. Some data are lost because 
certain antpols have extremely poor system tem- 



peratures, for various reasons: failed cryogenics 
systems, damaged LNAs, or offset pointing due 
to imbalanced electrical connections. More robust 
replacement hardware will make these issues much 
rarer. 

3.2. Calibration 

The calibration of our data is challenging be- 
cause of the modest u-v coverage, short integra- 
tion times, and wide temporal and spatial separa- 
tion between snapshots entailed by our observing 
mode. We will provide an overview of our cal- 
ibration process with more detailed descriptions 
presented in the subsections below. 

The snapshot observations of the calibrator 
sources were first self-calibrated to point-source 
models. All observations of the same calibrator 
and frequency were then imaged jointly to form a 
self-calibrated reference image with u-v coverage 
and sensitivity superior to those obtained in any 
one snapshot. The individual calibrator snapshots 
were then recalibrated to their corresponding ref- 
erence images. 

The snapshot observations of the science sources 
then had appropriate calibration solutions applied 
to them. Some snapshots (of both science targets 
and calibrators) were duplicated and had several 
different calibration solutions applied to allow for 
cross-checks of the calibrator flux models. Ref- 
erence images were then generated from the cali- 
brated observations and once again the individual 
snapshots were recalibrated to their corresponding 
reference images. 

Given the nature of this work, accurate broad- 
band calibrator flux density models are vital. The 
1999.2 analytic models listed in the VLA Cali- 
brator Manual were initially adopted and found 
to yield good results. These models were, how- 
ever, almost a decade old at the time of this work, 
and long-term VLA monitoring indicates that the 
fluxes of the calibrators evolve slightly on this 
timescale. Section lS^ describes consistency checks 
used to verify and fine-tune the flux density scale 
used in this work. 

3.2.1. Self- calibration 

The calibration process started with self- 
calibration of the calibrator observations. Band- 
passes and complex antenna gains in the indi- 
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vidual snapshots were first determined via the 
MIRIAD task smsmif cal, which assumes a point- 
source model. Spectral smoothing was used to 
increase the signal-to-noise ratio (SNR) in the 
high-frequency observations. 

System temperature (TLys) information was 
then computed and inserted into the datasets. 
This step was necessary because the ATA did 
not have online system temperature calibration 
at the time of this work. A Python script called 
calctsys was developed to estimate Tgys values 
from the variances in the real and imaginary parts 
of the visibilities across the band. The variances 
were used to compute per-basepol Tgys values as- 
suming an antenna gain of 153 Jy/K. Per-antpol 
Tsys values were then computed by performing a 
least-squares minimization assuming that the Tgyg 
of each basepol was the geometric mean of the 
Tsys of its two contributing antpols. The values 
obtained in this way had good reproducibility over 
the course of a segment and agreed well with lunar 
measurements made by ATA staff (R. Forster, pri- 
vate communication). The system temperatures 
could vary strongly as a function of frequency, 
with most, but not all, antpols demonstrating su- 
perior performance at lower frequencies. 

Reference images were then created from the 
group of calibrator snapshot observations. All of 
the snapshots of each source and frequency were 
merged into a temporary dataset that was imaged 
using natural weight ing. Multifrequency synthesis 
( Conwav et a/.|[l"990 l) was used to avoid bandwidth 
smearing and provide better u-v sampling. The 
images were deconvolved with the MI RIAD im- 
pleni entation of the CLEAN algorithm ([Hogbom 
1974I ). which incorp orates t he alg orith mic adjust- 
ments suggested by IClarkI (|l98(i ) and ISteer et al\ 
(jl984l ). After an initial round of imaging, the 



merged visibilities were self-calibrated to the gen- 
erated CLEAN component model and imaged once 
more. The second self-calibration was phase-only 
if the SNR of the data was less than 0.3. Testing 
showed that further iterations of calibration and 
imaging did not lead to improvements in the im- 
age quality. Figure |4] shows the reference image of 
NGC 253 at 1.3 GHz and an image generated from 
a single snapshot observation of the same source 
at the same frequency. Figure [5] shows reference 
images of M82 at four frequencies. The individual 
calibrator snapshots were then recalibrated using 



the CLEAN component model of the reference im- 
ages. 

3.2.2. Secondary calibration 

The complex antenna gains and bandpass so- 
lutions were copied from each calibrator snapshot 
to any observations made at the same frequency 
within 40 minutes of the calibrator observation, 
including observations of different calibrators. To 
allow further consistency checks, if multiple cal- 
ibrator observations could be applied to a given 
snapshot, it was duplicated and a different cali- 
bration was applied to each copy. The copy with 
the smallest time differential between the two ob- 
servations was treated as the "ofhcial" version of 
the observation, with the other copies not being 
included in the final analysis. The median time 
differential of the "ofhcial" observations was 21 
minutes, with values ranging from 3 to 38 minutes. 
The stability of the ATA over such timescales was 
qualitatively very good. Such stability is particu- 
larly noteworthy because of the large angles on the 
sky separating the calibrator and target sources 
(see, e.g., Table [2]). 

Once the observations of the science targets 
were flux and bandpass calibrated, system tem- 
peratures were computed and inserted into the 
datasets using calctsys. If the science targets 
had bright, off-center emission, the calctsys tech- 
nique would have computed incorrect Tgys values 
due to significant visibility variations across the 
passband, but this was not an issue for our se- 
lected sources. 

Reference images of the science targets were 
then generated from the calibrated observations 
in a manner similar to that described above. A 
reference image was generated for each observed 
combination of source, sky frequency, and calibra- 
tor. Before imaging, the merged data were phase 
self-calibrated against a primary-beam-corrected 
model gener ated from the NRAO VLA Sky Sur- 
vey (NVSS; ICondon et al] (|l998f)) catalog. The 
purpose of this step was to compensate for oc- 
casional antenna-based phase glitches that could 
occur between the calibrator observation and the 
observation it was applied to. To roughly cor- 
rect for the fact that the NVSS data only com- 
prise a single frequency, each NVSS model im- 
age was scaled by Sy{v)/ Sy{l.A GHz) before the 
phase self-calibration, where v was sky frequency 



7 



of the ATA observations and Su{v) was a simplis- 
tic model of the spectrum of the source being ob- 
served. The NVSS catalog sources corresponding 
to NGC 253 was found to poorly reflect its ex- 
tended structure, so a bootstrapped model of that 
source's spatial emission was used rather than the 
NVSS model. After this first round of phase self- 
calibration, the u-v data were imaged and recali- 
brated as described above. 

3.3. u-v modeling 

Fluxes were determined by fitting the u-v data 
with point-source and elliptical Gaussian mod- 
els based on the reference images. Our measure- 
ments are listed in TablesISHS The calibrators and 
Arp 220 were all well-modeled as point sources. 
M82 was modeled as an elliptical Gaussian. The 
deconvolved source sizes determined for M82 var- 
ied slightly over the range of frequencies observed. 
The major (minor) axis of M82 measured approxi- 
mately 60 (40) arcsec at 1 GHz and 35 (20) arcsec 
at 7 GHz, with a roughly constant position angle 
(PA) of ~80° (east from north). The higher- 
frequency values are consistent wit h observations 
made at 1.48, 4.87 32. and 92 GHz dSeaquist et al. 



Klein et al. 1988t Carlstrom fc Kronberg 



1985 

19911 ). The lower-frequency values are likely bi- 
ased by the large (~250") synthesized beam of the 
ATA at low frequencies. 

NGC 253 was modeled as the superposition 
of an elliptical Gaussian and a more compact 
source. At lower frequencies, the field of view con- 
tained bright sources away from the pointing cen- 
ter which were also modeled. The size of the ex- 
tended emission is strongly frequency-dependent, 
having a major (minor) axis of 480 (200) arcsec at 
1 GHz and 200 (50) arcsec at 7 GHz, with roughly 
constant PA of ~50 ° . The morphology of this re- 
gion at 1.5 G Hz is consistent with the "inner disk" 



dentified in Hummel et al. One should 



bear in mind that the spectra presented in Ta- 
bleland FigurcfTOl manifest not only the changing 
brightness of this region but its changing size as 
well. The more compact source was unresolved at 
lower frequencies and modeled as a point source. 
Above 3.5 GHz, however, it was possible for some 
structure in the compact source to be discerned, 
and it was found that an elliptical Gaussian with 
major (minor) axes of ^--^16 (~10) arcsec at a PA 
of ~60 ° was a practical model. This morphology 



is consiste nt with the "n u clear peak" seen in the 
near IR bvl Scoville et al\ (|l985h . 

The residuals to the fitted u-v models were im- 
aged to check for missed sources, poorly calibrated 
baselines, and so on. Further flagging was usually 
sufficient to solve any issues. Some datasets with 
severe RFI were not recoverable and had to be 
discarded. 

Figure[5]shows a histogram of the fractional flux 
residuals of the calibrator snapshots as compared 
to the analytic calibrator models. The calibra- 
tor fluxes are those derived from u-v modeling of 
calibrator snapshots that have been self-calibrated 
to the calibrator reference images. The vast ma- 
jority of the measurements agree to within 1%, 
demonstrating that the reference-image calibra- 
tion pipeline propagates fluxes with high fidelity. 

3.4. Uncertainties 

The MIRIAD u-v fitting task, uvfit, reports 
statistical uncertainties, but these fail to account 
for systematics and are only lower limits. We aug- 
mented these uncertainties with both relative and 
absolute terms; 

= sj + {i^S^^.f + {<Ja-l Jy)', 

where (Ji is the uncertainty used for the i**^ mea- 
surement. Si is the uncertainty for that measure- 
ment reported by the fitting routine, S^^i is the 
value of the measurement, and rj and a a are tun- 
able parameters. 

The choice of these parameters was derived 
from the statistics of the ensemble of observa- 
tions. For each observed combination of source, 
frequency, and calibrator, the weighted mean of 
the contributing observations was computed: 



Let there be N observed combinations of source, 
frequency, and calibrator, with rij observations of 
the J*'' combination. Let Sv,j be the weighted 
mean of the rij observations, and Si,.j^i and aj^i 
the values and uncertainties of the individual ob- 
servations. If the uncertainties have been correctly 
assessed and the measurements are drawn from 
normal distributions, the sum of the squares of 
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the normalized residuals, 



N rij 

EE 

j=i 1=1 



should be taken from a xt distribution, where k = 
SjLi "-i — ^ is the number of degrees of freedom. 
The expectation value of such a distribution is k. 

The value of was computed on a grid of 
r] and CTq for three groups of measurements: the 
full ensemble, those with < 1 Jy ("faint"), and 
those with 5*^ > 3 Jy ( "bright" ) . Figure [7] shows 
contours of = k (i.e., uncertainties agreeing 
with observed scatter) in the rj-aa plane. As ex- 
pected, the contours of constant are primar- 
ily along the rj axis for the faint group (for which 
CTa dominates the additional uncertainty) and pri- 
marily along the aa axis for the bright group (for 
which r] dominates the additional uncertainty). 
The = k contours for the three groups inter- 
sect at an approximate location of ?7 = 0.030 and 
(Ta = 0.014, i.e. an additional 3.0% fractional un- 
certainty and 14 mjy absolute uncertainty. These 
are the values of 77 and CTq that were adopted. The 
value of (Tq is approximately equal to the thermal 
noise in good conditions and is about half of the 
thermal noise at high frequencies, where Tgys gen- 
erally increases. 

If the uncertainties are correctly assessed and 
measurements are drawn from normal distribu- 
tions, the sum of the normalized residuals of each 
measurement is drawn from a normal distribution 
with mean and variance = '^jLi n-j. For the 
chosen parameter values, the sum of the residu- 
als was 0.52(T. This moderate value suggests that 
the adopted uncertainties characterize well the ob- 
served scatter in the measurements and that the 
weighted-mean approach does not mask large vari- 
ations in the individual measurements. 

3.5. Spectra 

We constructed broadband, high-resolution 
spectra out of the flux measurements described 
in ^13. 31 Repeated measurements were combined 
by taking a weighted mean as described in the 
previous section. Twenty-five percent of the spec- 
tral data are based on a single observation, while 
63% are based on at least three observations. The 
maximal number of contributing observations is 



13, of M82 at 5.0 GHz. The uncertainty of each 
spectral point is the propagated uncertainty of the 
weighted mean. 



The spectra are shown in Figures ISlfTT] along with 
measurements from the literature and model re- 
sults as discussed in ^ 

3.6. Self-consistency Checks 

We performed no absolute flux calibration of 
the ATA. The measurements we have made are 
thus measurements of flux density ratios, with ab- 
solute flux densities being derived from a priori 
spectral models of the calibrator sources. We can- 
not verify the models themselves, but we can as- 
sess whether our measured flux ratios are consis- 
tent with them. 

Let the "true" flux density of source X (i.e., the 
value that would be measured by a perfect, abso- 
lutely calibrated antenna at the time of observa- 
tion) be Tx , the a priori modeled flux density of X 
be Mx, and the flux density of X observed by the 
ATA after calibration to source Y be Ox,y ■ AH of 
these values vary with frequency. In the simplest 
possible model, at a flxed frequency Ox,y and Tx 
are directly related by a calibrator-dependent gain 
factor gy'. 

Ox,Y — gyTx- 

The gain is self-calibrated such that Oy.y = My, 
so that ideally 

Ox,y ^ My^ = MyRj,y, 
ly 

where y the true flux density ratio of X 
to Y. Not knowing whether this simple model is 
correct, we will define the observed flux ratio as 



Ox,Y 

My 



If this simple calibration model holds, the ob- 
served ratios should demonstrate consistency and 
closure: 



R 



X,Y 



■R 



Y,X 



= 1, 



R 



X.Y ■ ^Y.Z ■ 



R^ 



z.x 



(1) 

(2) 
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Furthermore, if the spectral models are accurate, 
it should be found that 



nO _ nM 



Mx/My 



(3) 



(If Mx is not known precisely, two measure- 
ments Y ^-iid Rx z still be used to check 
the model ratio Ryz by computing Ry z — 
Rxz/^XY-) If three of these properties are 
observed, there is good reason to believe that 
RxY = RxY ^^'^ Ox.Y = Tx, though it is im- 
possible to demonstrate the latter without abso- 
lute flux measurements. 

We compute these ratios and plot them in Fig- 
ure 1121 Our measurements do not have enough 
precision for us to believably solve for spectral 
models of the calibrators from scratch. We do, 
however, solve for a change in the overall scaling 
of the VLA 1999.2 models, i.e., for three parame- 
ters /48, /i47, and /286 such that Mx = fxVx, 
where Vx is the VLA model of source X and 
we denote the calibrators by their 3C numbers. 
Given the constraints implied by Equations [T] and 
[21 there are only two independent Rx y, so a nor- 
malizing assumption must be applied to derive 
the three fx from observations. We choose to 
fix /48/147/286 = I7 which is equivalent to requir- 
ing that the total relative change to the models. 



In^ /48 



In^ /i47 -|- In^ /286J be minimized. 



To solve for the fx from the data, we define 
a = /48//147 and b = /48//286- Minimizing the 
logarithms of R'^ y/^x y in a least-squares sense 
leads to a = 1.018 and b — 1.007 with a root 
reduced parameter {\/^) of 1.07. This mini- 
mization included two precise measurements made 
at the VLA in 2008 September: R%.2S6 = 1-066 ± 
0.003 at 1.465 GHz and i?^8.286 = 0.7349 ± 0.002 
at 4.885 GHz (R. Perley, private communication). 
Using our normalizing assumption, we find /48 — 
1.011, /i47 0.989, and /286 = 1-001. This is 
consistent with evidence from the VLA that the 
flux density of 3C 48 has recently risen slightly 
(R. Perley, private communication). 

Our measured flux ratios are well-described 
within their uncertainties by the rescaled VLA 
models with the exception of the i?i47 286 data. 
These observations were made in June 2009, well 
after the main observing for this work was com- 
pleted, after it was realized that there were only 
a few data comparing these calibrators. Unlike 



the 2008 observations, these occurred in the day- 
time, and solar interference was a signiflcant prob- 
lem with the low-frequency data. Nevertheless the 
high-frequency ratios, for which solar interference 
was much less relevant, disagree significantly with 
the adopted models. We have not been able to 
identify the cause of this disagreement, but we do 
not believe that it is due either to genuine inaccu- 
racy in the models or to problems in the calibra- 
tion and u-v modeling pipeline. If the former were 
the case, it would represent an unprecedented vari- 
ation in the flux of some combination of 3C 286 
and 3C 147. If the latter were the case, the con- 
sistency property of Equation [T] would not hold as 
well as it does in the observations. We believe the 
robustness of the rest of the results provide ample 
reason to believe that the models and reduction 
pipeline are generally robust and accurate. 

4. Spectral Modeling 

In order to assess the steepening of the non- 
thermal components of the observed spectra, we 
modeled them wit h a two-component model simi- 



lar to that used in lNiklas et all ( 199711 : 



(4) 



where S'th(i^) is the thermal contribution and 
Sntii^) is the nonthermal contribution. Free- free 
absorption (FFA) is likely to play a signiflcant role 
in determining the spectra of our target galaxies. 
We attempt to account for its effects by including 
as a parameter in our model ti, the free- free op- 
tical depth at 1 GHz. The thermal component is 
then 



5thH = S,(Te)r!(l-e-"«), 



(5) 



where Tg ~ lO'' K is the temperature of the ther- 
mal electrons, is the angular size of the emitting 
region, and rg is the free-free optical depth as a 
function of frequency. We use the canonical fre- 
quency dependence of the free-free optical depth, 



-2.1 



(6) 



where i//GHz — Vg. In typical cases, the large- 
scale thermal emission in gal axies at cm wave- 
lengths is optically thin [e.g. iDuric et~ai. 19881: 
Price fc Duridll992l ). rg ^ 1, in which case the 



thermal emission has the well-known frequency de- 
pendence of 5'th oc v~'^'^ . The galaxies we consider 
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are, however, decidedly atypical, with dense, com- 
pact emission regions, so we avoid this approxima- 
tion in our modeling. Our results indicate, how- 
ever, that it broadly holds even for the galaxies in 
our sample (cf. Figures [51 fTT|) . 

For simplicity we assume that the thermal and 
nonthermal components are well-mixed and hence 



If this i s not the case (as is 
Beck fc Reichlll985l 



share the same 
true of the Milky Way, e.g 
this assumption would require loosening in a more 
detailed analysis, although at small values of 
its impact is minimal. In the well-mixed case, the 
nonthermal emission is attenuated so that 



1 



Tff 



(7) 



where ^nt is the unabsorbed nonthermal spectrum. 
Because our fits model the thermal and nonther- 
mal contributions simultaneously and we do not 
approximate any FFA-related effects, we should 
obtain good results even when ti is large. We 
model the unabsorbed nonthermal component as 

5'ntK) = 10'4+Slogio''.<,+C[logio^.,l'jy^ (8) 

a parabola in log-log space. Our parabolic spectral 
model is not intended to be well-motivated phys- 
ically, but rather to be simpl e and fl exible. The 
theoretical spectra modeled in lT-l-06l to which we 
compare our results in 53 are derived numerically 
and hence do not suggest a more approp riate an- 
alytic expression for us to employ. While ICondonI 
(|l992l ) and others typically use broken power laws 
to model aged synchr otron spectra, one of the key 
contentions of is that other loss processes 

will signiflcantly modify raw synchrotron spectra 
in starbursts. We note that functional forms anal- 
ogous to ours have proven t o be good e mpiri- 



cal models in works such as iBaars et all (|l977l ) 



and the VLA broadband calibrator models. The 
parabolic model has the additional benefits of be- 
ing linear in its parameters (in log-space) and be- 
ing reducible to a pure power-law spectrum (ig- 
noring FFA) simply by fixing C = 0. 

An alternative model for the effect of FFA upon 
the nonthermal emission is one in which the ab- 
sorption screens but is not mixed with the non- 
thermal emitting medium. 



We disfavor this model for our targets because 
they are starbursts, and their diffuse synchrotron 
emission propagates through a compact ISM rid- 
dled with H II regions and supernova remnants 



McDonald et aZ.ll2002HRodriguez-Rico et al. 



2004 . regarding M82) . In the analysis that follows 
we will, however, occasionally compare results for 
the two absorption models. 

4.1. Isolating the thermal emission 

Excepting those with active nuclei, most galax- 
ies have a nonthermal spectrum much steeper than 
the aforementioned a ~ —0.1 of optically-thin 
thermal free-free emission. Because the nonther- 
mal emission dominates at cm wavelengths, good 
higher-frequency data are needed to usefully con- 
strain the thermal component. However, above 
v ~ 100 GHz, thermal emission from cold dust be- 
gins contributing non-negligibly to the total fiux. 
We incorporated fiux densities from the literature 
with 20 < I// GHz < 100 into our fits for this rea- 
son. Table [6] summarizes the parameters of some 
of the relevant observations found in the litera- 
ture. Those observations, and others not listed in 
table [6l are shown in Figures [8lfTT] for comparison 
to our observations. 

Comparing fluxes in this way requires care 
because interferometers resolve out spatially ex- 
tended emission. The largest angular scale (LAS) 
to which an interferometer observation is sensi- 
tive is inversely proportional to its shortest base- 
line as measured in wavelengths. In a very real 
way, two different interferometer observations, or 
an interferometer observation and a single-dish ob- 
servation, are not "looking at" the same source. 
We deal with these considerations on a source-by- 
source basis below. 

4.1.1. M82 

M82 is well-observed in the 20-100 GHz range 
with both interferometers and single dishes. The 
interferometric fluxes are systematically smaller 
than the single-dish fluxes by about 25%, with the 
exception of the unus ually low 87.2 GH z single- 
dish measurement of Jura et 

'd\ (|l978f) . Unlike 



(9) 



the other single-dish data, which come from maps, 
this observation was made with a single pointing of 
an 11 m dish (HPBW 75"). It furthe rmore needs 
conversion to the flux density scale of lBaars et al. 
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(jl977f ). We regard it as an outlier and discard it. 
All of the interferometer observations we con- 



sider |Carlstrom fc Kronberdl991: Ro dn'guez-Rico et al 
2004 ICarlstrom fc KronberdllQQll: ISeaquist et ail 



19961 ) have LASs of ~50". The deconvolved size 
of the radio emission of M82 is ~35"x -10" 
from cm wavelengths (this work) up to 92 GH z 



(|Klein et all Il988t ICarlstrom fc Kronberd Il99ll ). 



so significant flux from M82 should not be resolved 
out; the source of the excess emission detected 
by the single-dish observations is unclear. For 
safety, we compare the ATA fluxes only to other 
interferometric measurements. This criterion is 
coincidentally equivalent to the selection of only 
post-1990 observations. 

4.1.2. NGC 253 

The extended spatial structure of the emission 
of NGC 253, as detected by the ATA, makes the 
modeling of its spectrum more difficult. We do 
not expect the high-frequency interferometer ob- 
servations of NGC 253 listed in Table [6] to be 
sensitive to the extended component because of 
its steep spectral index and the constrained LASs 
of the observations. Simple simulations of high- 
frequency VLA observations of the galaxy sup- 
port this conclusi on. Two of the FIR observations 
of NG C 253 (Rieke et fl/.lll973l : iHildebrand eTd. 



19771 ) are made with single dishes with ~65" 
beams that could be sen sitive to th e exten ded 
emission region; however, IWeifi et al. (l2008h re- 
solve the central component from the extended 
source and find that emission from the former 
dominates the latter. Because of these facts, our 
spectra link the FIR and high-frequency measure- 
ments with fluxes for only the central component 
of NGC 253. 

4.1.3. Arp 220 



Th e observations o f Arp 220 made byllmanishi et al. 
(|2007l ). lNorri^ (| 19881 ) and lZhao etld] (|l996h have 
LASs that are too small for our purposes and 
are not used in o ur fits. The 97 GHz single- 
dish observation of Anantharamaiah et al. ( 2000l ) 
is significantly discrepant from 87 GHz observa- 
tions made with the Plateau de Bure Interferome- 



ter (jRadford et a/.lll991l). the N obeyama Millime- 
ter Array ( Imanishi et al. 2007t ). and a 110 GHz 
measurement made with the Owens Valley Radio 



Obse rvatory (OVRO) mm Array (jScoville et al. 
199l[ ). The 113 GHz measurement also pre- 
sented in lAnantharamaiah et al ] (|2000[ ) clusters 
much better with these other observations and 
a power-law extrapolation of FIR observations 
of A r p 220 dCarico et al\ [l992: D unne & Ealei 



dea 

m 



20011: lEales et al\ 119891 : iRigopoulou et al\ Il99 
Woodv et al. fEiii). Once again, we use the in- 
terferometric observation but not the single-dish 



4.2. Modeling results 

We performed least-squares fits between the av- 
eraged ATA data and the chosen high-frequency 
measurements for each source, using a Levenberg- 
Marquardt algorithm with constrained to be 
nonnegative. In the case of M82 we used = 2400 
arcsec^, its modeled angular size at 1 GHz. For the 
core component of NGC 253, we used its mod- 
eled size as resolved at higher frequencies, 160 
arcsec^. The ATA d id not resolve Arp 220 but in 



Scoville et al\ (|l99ll ) it is reported that the emis- 
sion comes from an extended component of ^1" 
X ^^15" and a stronger central component of —2 
arcsec^. We adopt an intermediate value of 17 
arcsec^. In all cases we used Tg — 10^ K. The 
numerical results of our fits are reported in Ta- 
ble [71 The uncertainties shown are marginalized 
values determined via Monte Carlo simulations 
with 10000 realizations for each source. To pro- 
vide a more physical interpretation of B and C, 
we also transform them, respectively, into 05, the 
spectral index of the nonthermal model compo- 
nent at 5 GHz, and A = Aa/dilogi^v) at the 
same frequency. We also compute /th, the ther- 
mal fraction (not correcting for FFA) at 1 GHz. 
We computed 3a confidence limits in the total flux 
of the best-fit models at 40 frequencies logarithmi- 
cally spaced between 0.7 and 100 GHz. Finally, we 
computed Icr confidence regions in the B-C plane 
of parameter space. This was done by tracing the 
contour of = Xmin + 2.3, with xLin being the 
minimum achieved in the overa ll fit a nd the 
offset parameter coming from lAvnil (|l976l ). The 
value of assigned to each location in the B-C 
plane is the minimum value of attained in a fit 
with B and C fixed and A and ti free. 

The parameter values that are obtained de- 
pend strongly on the highest-frequency data 
points. Fits that were run without any very 
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high-frequency {ly > 50 GHz) data points often 
converged to the physically-suspect lower bound 
of Ti = 0. To investigate the possible effects of 
inaccurate high-frequency observations on our re- 
sults, we ran a set of Monte Carlo simulations 
in which the uncertainty associated with each 
high-frequency measurement was doubled. In this 
case, the resulting uncertainties in the modeled 
parameter values increased by a factor of ~1.5, 
demonstrating the extent to which the results are 
affected by those few data points. While we have 
done our best to choose our high-frequency data 
points carefully, it remains true that our results 
are dependent on the choice of which measure- 
ments to include and, of course, the accuracy of 
those measurements. 

4.2.1. M82 

The M82 fit has y/^ ^ 1.21. The ATA data 
and fit results are presented in Figure[8l The value 
of Or that we obtain agrees well with that found 



(|l98ah . a = -0.68 ±0.02. We, how- 
ever, find a higher thermal content at 32 GHz, be- 
cause the possibility of a steepening nonthermal 
spectrum allows for a larger thermal contribution. 

A fit with a fixed value of C = (i.e., a pure 
power-law nonthermal spectrum) had a compara- 
ble v^xf = 1.34, but also ri = 0, a physically 
implausible result. In the C ~ case, the many 
ATA data points tightly constrain the power-law 
slope of the nonthermal component, but the high- 
frequency points lie significantly below the non- 
thermal line. The thermal component, being spec- 
trally shallower than the nonthermal component, 
can only make this situation worse, so the fit con- 
verges to a solution with no thermal component at 
all. If any thermal component is present, the non- 
thermal spectrum must steepen so that the high- 
frequency points can be modeled successfully. 

4.2.2. NGC 253 

The NGC 253 fit has = 1-02. The ATA 
data and model results for the core component of 
NGC 253 are presented in Figure [5] The sym- 
bols used are analogous to those in Figure [H Fig- 
ure [in] shows the extended and total fluxes in con- 
text with total-flux data from the literature. The 
same models as in Figure IHl are overlaid to ease 
comparison. 



A fit with a fixed value of C = had \/x^ = 
1.51 and ti = 0.25. Here, a change in the spec- 
tral index is detectable in the ATA data alone, so 
when the nonthermal component has a constant 
slope, Ti must increase to induce more curvature 
at low frequencies via FFA. In this case, once again 
the high-frequency data points lie well below the 
points attained by the model. 



4.2.3. Arp 220 



Xr = 1-23 and ri = 



The Arp 220 fit has 
and is shown in Figure fTTl This result is due to the 
dramatic steepening of the spectrum of Arp 220 at 
ly ^ 10 GHz. A purely nonthermal fit that accom- 
modates the high-frequency observations must un- 
derestimate the low-frequency (i/ < 2 GHz) data 
points significantly. The addition of any thermal 
contribution would require an even steeper non- 
thermal component to match the high-frequency 
measurements, which would make the discrepan- 
cies on the low end even worse. Adjustments to 
the assumed angular size of Arp 220, which affect 
the overall scaling of the thermal component, do 
not improve the outcome. We report the results 
of our fit in Table [7] but caution that they are un- 
satisfactory. We exclude the Arp 220 results from 
further analysis in the following sections. 

If the nonthermal curvature parameter C is 
fixed to zero, the model cannot accommodate the 
data, with y/)(^ — 2.77. This result quantita - 
tively confirms the supposition of I Condon ( 19921 ) 
that the nonthermal spectrum of Arp 220 steepens 
in the cm regime. Our data cause serious prob- 
lems for the empirical spectral model of Arp 220 
presented in lSopp fc Alexanden (|l99ll ). which pre- 
dicts a much steeper spectrum at frequencies of a 
few GHz than we find. Their model was incon- 
sistent with several measurements at ~ 2 GHz 
shown in the paper, and our ATA data make it 
clear that those measurements are generally accu- 
rate and the simple power-law used in that work 
is insufficient to match the data. 



5. Discussion 

We test whether our observational results for 
M82 an d NGC 253 are consistent with the predic- 
tions of iT-HOel by comparing the curvature of the 
FFA-corrected nonthermal spectra with the out- 
put of the numerical model presented in that work. 
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We obtained an updated version of the relevant 
code (E. Quataert, private communication) and 
modified it to allow more convenient tuning of the 
input parameters. These parameters are (1) Eg, 
the surface density of the gas disk; (2) /left, the 
effective scale height of the gas disk; (3) p, the 
power-law index of the injected CR energy spec- 
trum, dn (X E~PdE; and (4) 77, the ratio of the en- 
ergy density of the magnetic field to the pressure of 
the ISM, = Ty-STT^GS^. The value of hes affects 
the characteristic timescales of CR cooling via ion- 
ization and free-free losses by setting the density 
of the neutral ISM as encoun tered by the CRs, 
ncR = Sg/2mp/ieff- (In IT-I-OgI hcs is specified by 
two parameters, the physical scale height of the 
gas disk h and an ISM/CR interaction efficiency 
parameter /: hcs = h/ f. They are degenerate for 
our purposes.) Because the mean value of the ISM 
density may not correspond well to the value of 
ncR due to clumping and other effects, hes may 
likewise not correspond well to the actual scale 
height of the gas disk. We take the surface density 
of the disk of M82 to be 0.081 g cm"^ and that of 
NGC 253 to be 0.067 g cm'^ (jThompson et al\ 
I2OO9L and references therein). Note that these 
recently-derived value s are signif icantly different 
than the ones given in iKennicuttI ([l998). 



ISM pressure are in equipartition (jParker Il966h . 



Given a choice of S^, hcff, p, and 77, we used a 
simple least-squares fit to determine values for the 
spectral parameters B and C. We used the 
model to compute flux densities at 25 frequencies 
logarithmically spaced between 1 and 7 GHz and 
then fit a bent power law (as in Equation[5]) to the 
data. All of the spectra generated in this way were 
very well-described by such a model. The derived 
parameters were not sensitive to the details of the 
choice of frequencies at which to sample the model. 

We computed such B and C values for the 
T-I-Od model on two grids of parameters. The first 



had Sg set to the observed M82 value, h^s ranging 
between 10 and 250 pc, p ranging between 2 and 
2.5, and 77 ranging between 0.02 and 1, all step- 
ping through 10 values logarithmically. The lower 
bound of 77 is the approximate value for which the 
synchrotron cooling time, Tgyn, and the timescale 
for escape from the galaxy via advection by a 
galactic wind, Tosc, are equal. For 77 (and hence 
B) smaller than this, the assumptions of the lr+Od 
code are no longer upheld. The largest plausible 
J] is ~1, at which point the magnetic field and 



The second grid was identical except that Eg was 
set to the NGC 253 value and the lower bound 
of 77 was 0.04, the analogous limiting value. Fig- 
ure[T3]shows the values of B and C attained by the 
model and compares them to the la and 2a con- 
fide nce inte rvals associated with our observations. 
The lT-H06l models can accommodate a wide range 
of spectral indices and exhibit spectral steepening 
for all parameter values. Our observed spectra, 
however, exhibit more steepening (more negative 
C) than the models. The models are consistent 
with observations of both galaxies at the 2a level 
but only with M82 at the la level. If the source 
of the FFA is configured as a screen in front of the 
nonthermal emission rather than a medium mixed 
within it (Equation [9]), the C values derived for 
our observ ations are less negative, ~ —0.10, and 
the IT-I-06I models have large areas of overlap with 
the observations of both galaxies even at the la 
level. The values of the fits in this case are 
virtually identical to those of the fits in the well- 
mixed case. 

We quantify this agreement in Figures [14] and 
I15[ which show the r egions in the rj-hcff plane 
where the results of the T -1-061 model are consistent 
with the observed properties of M82 and NGC 253, 
respectively, at the 2a level, for p = 2 a nd 2.3. 
Because our data, unlike the outputs of the IT-I-06I 
models, are affected by FFA and thermal emis- 
sion, they cannot be compared to the models di- 
rectly without taking into account these effects. 
To do so requires some sort of spectral modeling, 
so we compare the spectral curvature parameters 
B and C found for o ur dat a with ones derived 
for the outputs of the T -1-061 models. (If we cor- 
rected our observations for the modeled FFA and 
thermal emission and then compared them to the 
T-^06l models, we would be effectively hiding the 



model-dependence of the correction and discard- 
ing the information that the data have been cor- 
rected assuming an underlying parabolic nonther- 
mal spectral model.) We assumed a bivariate nor- 
mal probability distribution for the observed val- 
ues of B and C, with the parameters of the distri- 
bution estimated by applying the maximum like- 
lihood method to the results of the Monte Carlo 
simulations. The agreement between the shape of 
the analytic distribution function and contours of 
constant Ax^ was very good. Using the bounds on 
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hes and rj given above, p = 2.5 is rejected at the 
99.96% confidence limit for M82 and the 99.99% 
confidence limit for NGC 253; larger values of p 
become acceptable if the lower limit on heS is de- 
creased. For NGC 253, p = 2.3 is also rejected at 
the 99.99% confidence limit. 

For M82, virtually all physically plausible val- 
ues of 1] are allowed. For p — 2, the allowed val- 
ues of hcff are of order 100 pc, broadly consistent 
with observed radio sizes. For p — 2.3, agree- 
ment requires hcs ~ 20 pc. Such sizes have been 
proposed for the neutral gas disks o f some ultralu- 



the p — 2.3 case. A precise independent measure- 
ment of p cou ld thus provide a valuable test of this 
aspect of the model. 



minous infrared galaxies (ULlRGs; IScoville et al. 



1997: Carilh et al 



1991), but are unlikely to be 
found in less extreme galaxies. Furthermore, if 
the CRs interact with the ISM at somewhat less 
than it s nominal mean density (the parameter / of 
T-HOel is <1), as is sugg ested by meas urements of 
CRs in the Milky Way (|Connelllll998h . the actual 
scale height of the disk becomes even less phys- 
ically plausible. As rj increases, the allowed h^s 
decreases for both values of p. This is because the 
stronger magnetic field steepens the synchrotron 
spectrum, requiring a smaller h^s to increase ncR 
and consequently the amount of spectral flattening 
due to ionization and free-free emission off of neu- 
tral hydrogen. As p increases, the allowed values 
of /left rapidly decrease for the same reason. 

The NGC 25 3 data are only marginally con- 
sistent with the iT-H06l model within the range of 
parameters we consider, and the statistically ac- 
ceptable results have hcs « 50 pc. If the source of 
FFA is configured as a screen rather than mixed 
with t he nonthermal emission, much more of the 
T-I-Od model parameter space is allowed, with a 
general configuration of allowed parameters simi- 
lar to that obtained for M82. 

The relationship between heS and the observed 
radio sizes of star burst galaxies /lobs is signific ant 
within the T -1-061 model. Estimates in T -1-061 in- 
dicate that if hcS ~ /lobs and B ^ -Bmin, free-free 
emission off of neutral ISM nuclei becomes a signif- 
icant sink of CR energy for > 1 g cm~^, which 
would disrupt the observed linearity of the FIR- 
radio correlation at high surface densities. Since 
the free-free timescale is independent of the mag- 
netic field while the synchrotron timescale goes as 
Tgyn oc B~^/'^, a stronger magnetic field, B ~ i3eq, 
avoids this inconsistency. Our observations have 
^cff ^ ^obs in the p = 2 case and hcs ^ /lobs in 



6. Conclusions 

Our observations of M82 are consistent with the 
scenario proposed in T-l-06 . Our observations of 
NGC 253, on the other hand, are only marginally 
consistent with it. We decisively detect the spec- 
tral steepening at ~ 1 GHz predicted in that 
work. The uncertainties in our spectral model- 
ing, due to uncertainties in the individual obser- 
vations, the difficulties of untangling the thermal 
and nonthermal emission, and the inaccuracies in- 
herent in our spatially-unresolved approach, make 
it difficult to si gnifican tly constrain all of the pa- 
rameters of the IT-H06I spectral model or to reject 
its entire parameter spa ce. We do, however, reject 
with high confidence the T -1-061 model with a steep 
electron energy injection spectrum power-law in- 
dex ofp = 2.5 for M82 and N GC 25 3. We also 
reject with high confidence the T -1-061 model with 
p = 2.3 for NGC 253. We are unable to place lim- 
its on the str ength of the magnetic field required 
by the lT-l-06i model with our current data. Fur- 
ther work to constrain the efficiency with which 
CRs interact with the ISM or the allowed values of 
p would significantly shrink the model parameter 
space. Observations of a larger sample of galaxies 
would allow a statistical treatment of the topics 
that we have dealt with here on a case-by-case ba- 
sis. It would likely also be fruitful to compare our 
detailed radio results with FIR observations. Such 
an approach would be a frequency-space analog 
to the recent works investigating the connections 
between FI R and radio em i ssion at high spatial 



20061) 



resolutions (Marsh fc Helou 1995; Paladino et al. 



The broadband, continuous spectra that we 
have obtained with the ATA are precise and re- 
peatable. The quality of our measurements is lim- 
ited by several factors. The u-v coverage of our 
snapshot observations is good but suboptimal for 
complex sources, which made the imaging of tar- 
gets such as NGC 253 more difficult. The sensitiv- 
ity the system is a limiting factor for snapshot ob- 
servations of fainter sources such as Arp 220. Both 
of these quantities were limited by the limited 
number of antennas available for observing during 
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the commissioning of the ATA-42. The planned 
expansion of the ATA to ~350 antennas will do 
much to improve both of them. The fact that the 
ATA was undergoing commissioning at the time of 
these observations meant that several important 
aspects of observing and data reduction, such as 
polarization calibration and RFI mitigation, were 
incomplete. Ongoing development since these ob- 
servations were made has led to improvements in 
virtually all aspects of the system performance. 
Of particular interest is the deployment of a sec- 
ond independently- tunable correlator, which dou- 
bles the already-fast throughput of our spectral 
observations and further improves the prospects 
for fast, precise, broadband spectral monitoring of 
sources. 
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Table 1 
Observations 



Date (UT) 


Segment # 


Sources 




Frequencies (GHz) 


20080912 


1 


3C 48, 3C 147, M82, NGC 253 


1.4, 


1.5, 2.0, 2.1, 3.1, 4.0^^ , 4.P , 








5 


0, 5.1 


20080913 


2 


3C 286, Arp 220 


1.5, 


1.6, 2.5, 2.6, 4.5, 4.6 




3 


3C 147, NGC 253 


1.1, 


1.2, 1.3, 2.6, 2.7, 3.6, 3.7^^ , 5.6, 




4 


3C 48, 3C 147, M82 


5 

3.2, 


7 

3.3, 4.4, 4.5, 5.5, 5.6, 6.6, 6.7 




5 


3C 286, Arp 220 


1.1, 


1.2, 2.5, 2.6, 5.0, 5.7, 6.3 


20080914 


6 


3C 48, Arp 220 


1.1, 


1.2, 2.5, 2.6, 5.0, 5.7, 6.3 




7 


3C 147, NGC 253 


2.7, 


3.1, 3.3, 4.5, 4.7, 6.1, 6.3 




8 


3C 48, 3C 147, M82 


1.0, 


1.1, 2.6, 2.9, 4.8, 5.0 




9 


3C 286, Arp 220 


1.9, 


2.9, 3.2, 3.3, 5.4 


20080924 


10 


3C 286, Arp 220, M82 


1.4, 


1.7, 1.8, 2.0, 3.4, 3.5, 3.6, 5.0, 








5 


2, 6.0 


20081015 


11 


3C 48, 3C 286, Arp 220 


0.5, 


0.7, 0.9, 1.4, 4.3, 4.4, 4.7, 5.0, 








6 


2, 6.7, 7.0 




12 


3C 48, 3C 147, M82 


1.4, 


2.4, 2.5, 4.3, 4.7, 5.0 




13 


3C 147, NGC 253 


1.4, 


1.7, 1.8, 1.9, 4.5, 4.6, 4.9, 5.0, 








5 


8, 5.9, 6.0, 6.4, 6.7, 7.0 




14 


3C 48, 3C 147, M82 


6.1, 


6.2, 6.3, 6.4 


20090307 


15 


3C 48, 3C 286 


1.4, 


3.1, 5.0, 6.0, 7.0 


20090621 


16 


3C 147, 3C 286 


1.0, 


1.4, 1.5, 5.0, 5.1, 5.5, 6.1, 6.5, 



7.0 



^Discarded due to severe RFI. 
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Table 2 

Snapshots made during segment #12 



Timp 


± tjiiitiiig, 




n ,1 i nn 

J— /IC V d L HJll 


7/ 1 


J- Ut/Uo 


CTTT'l 








\\J[±1.Zj J 




03:50 


3C 48 


77 


39 


1.4 


2.0 


03:53 




77 


40 


2.4 


3.0 


03:55 




78 


40 


2.5 




03:59 


M82 


359 


21 


1.4 


2.0 


04:02 




358 


21 


2.4 


3.0 


04:04 




358 


21 


2.5 




04:08 


3C 48'^ 


79 


43 


4.3 


5.0 


04:11 




80 


43 


4.7 




04:13 




80 


44 


5.0 




04:17 


M82 


359 


21 


4.3 




04:20 







21 


4.7 




04:22 







20 


5.0 




04:28 


3C 48'^ 


82 


47 


1.4 


2.0 


04:31 




82 


47 


2.4 


3.0 


04:34 




83 


48 


2.5 




04:38 


M82 


1 


21 


1.4 


2.0 


04:41 




2 


21 


2.4 


3.0 


04:43 




2 


21 


2.5 




04:47 


3C 147 


40 


19 


1.4 


2.0 


04:50 




40 


19 


2.4 


3.0 


04:52 




40 


20 


2.5 




04:56 


3C 48 


86 


52 


4.3 


5.0 


04:58 




86 


52 


4.7 




05:01 




87 


53 


5.0 




05:04 


M82 


4 


21 


4.3 




05:07 




4 


21 


4.7 




05:09 




4 


21 


5.0 





"^A blank entry indicates no change from the previous value. 

''Had it been above the 18° elevation limit of the ATA, 3C 147 would have been observed at these times. 
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Table 3 

FLT'X density MEARrREMENTS OF M82 



y (GHz) Flux (Jy) Calibrator 



1 n 

J- .(J 




"W, 48 




8 Qfi + fl 1 fi 


3C 147 


1.1 


c 70 + n o« 


'\C 48 






HC 147 


1.4 


7 oc; -1- n 1 1 


'\C, 48 




7 Q7 _(_ n in 


osfi 


1.5 


\J'\J\j —1— W . J- u 


3C 48 


1.7 


fi fin + fl OQ 


3C 28fi 


1.8 


fi Q4 _|_ n nn 


3C 286 


2.0 


fi Qfi _(_ n 1 Q 


3C 48 




fi 1 Q + n nQ 


"^C 28fi 


2.1 


c, Q7 _(_ n 1 Q 


■^C 48 

V-^ 4:0 




c: on _(_ n 1 


3r! 48 

40 


Zi.U 


(J .tJ(J HZ u. 1 1 


40 






'^C 147 




4 QQ _|_ n 1 1 


3r! 48 

k)\J 4;0 




4 QQ _|_ n 1 1 

'±.00 —1— L/.-L J- 


3C 147 

tJVj X4: J 


3.0 


4 QO 1 n nq 


3C 48 


3.1 


4 7n _|_ n 1 n 


3C 48 


3.2 


4 54 1 A n7 


3C 147 


3.3 


4 CA 1 n n7 


3C 147 


3 4 


4 44 _|_ n nfi 


osfi 


3 


4 QO _(_ A nfi 


"^C 28fi 


o.u 


^.±tj ZC L/.UU 






Q oc _L n nq 

O.OtJ ZC W.Ut? 


OVJ 4;0 


/I A 


0. ( □= U.UD 


op 1 /1 7 


4.5 


3.67 ±0.06 


3C 147 


4.7 


3.55 ± 0.08 


3C 48 


4.8 


3.53 ±0.08 


3C 48 




3.60 ±0.09 


3C 147 


5.0 


3.38 ±0.05 


3C 48 




3.50 ±0.07 


3C 147 




3.46 ±0.05 


3C 286 


5.1 


3.16 ±0.10 


3C 48 




3.52 ±0.11 


3C 147 


5.2 


3.34 ±0.05 


3C 286 


5.5 


3.34 ±0.05 


3C 147 


5.6 


3.20 ±0.06 


3C 147 


6.0 


2.96 ±0.05 


3C 286 


6.1 


2.99 ±0.10 


3C 147 


6.2 


2.91 ±0.10 


3C 147 



22 



Table 3 — Continued 



V (GHz) 


Flux (Jy) 


Calibrator 


6.3 


2.98 ±0.10 


3C 147 


6.4 


2.90 ±0.10 


3C 147 


6.6 


2.76 ±0.06 


3C 147 


6.7 


2.86 ±0.06 


3C 147 



23 



Table 4 

Flux density measurements of NGC 253 



V 




Flux (Jy) 




Calibrator 


(GHz) 


Core 


Extended 


Total 




1.1 




3 Q2 ± 07 


G qo + 08 


3C 147 


1.2 




Q 00 _(_ A AO 

0.00 —1— \J .\J<J 


c 7r, _|_ 1 


3C 147 


-L .O 


2 7fi + n 04 


Q Ki _(_ n A(? 


(? 29 + 07 


3C 147 


1.4 


_|_ n no 


Q 1 7 _|_ f) 1 n 


c: qq + 1 

tJ.t/tJ —1— U.J-O 


3C 48 






Q 22 _(_ A A7 


c; q7 _(_ 09 


3C 147 


1.5 




2 50 + n no 


c 1 1 1 Q 


3C 48 


1 7 

J. • \ 


Zi.*±0 —1— \J>\)\J 


2 fis + n 07 


5 1 fi + 09 


3C 147 




QQ _|_ n ric; 


2 c:o + 07 

z .<JU —1— yJ .\j t 


4 cq _(_ 08 


3C 147 


1 Q 

-L .t/ 


_|_ f) nc; 


2 42 + 07 


4 7c _(_ 08 


3C 147 




Qc: _|_ n nc 


2 2Q + 10 


4 G4 _|_ 12 


3C 48 


2.1 




1 Qfi + 09 


4 21 +012 


3C 48 


2.6 


X.C/tJ _1_ W«\JtJ 


1 00 -(- 05 


Q 75 _(_ 06 


3C 147 


9 7 


ni + n f)9 


1 74 + 03 


Q 7 c: _|_ 04 


3C 147 




1 Q« + n 07 


1 7fi + 09 


Q 74 _(_ 1 2 


3C 48 




1 «Q + n ofi 


1 70 + 08 


Q r;q _(_ 1 1 


3C 48 




1 s":! + n no 


1 CI _|_ 03 

J.. (J J. _l_ u.wo 


Q Q4 _l_ 04 


3C 147 


3.3 


1 7fi + n 02 


1 40 + 04 


1 fj 1 04 


3C 147 


3 fi 


1 7^; + n 04 


1 + 07 


Q 1 Q _(_ 08 


3C 147 


3 7 


1 71^ + n n» 


1 1 c _|_ 1 7 


2 qo + 19 


3C 147 


4 5 


1 _(_ A 02 


07 + 04 


2 47 + 04 


3C 147 


/I 1^ 

4.0 


1 /ifi -u n 0/1 

±.4J:D It U.U'i 


i.UD □= U.U ( 






4.7 


1.47 ±0.03 


1.00 ±0.09 


2.47 ±0.10 


3C 147 


4.9 


1.35 ±0.04 


1.01 ±0.07 


2.35 ± 0.08 


3C 147 


5.0 


1.37 ±0.10 


0.65 ±0.15 


2.01 ±0.18 


3C 48 




1.45 ± 0.04 


0.81 ±0.09 


2.26 ±0.10 


3C 147 


5.1 


1.34 ±0.07 


0.90 ±0.15 


2.23 ±0.17 


3C 48 


5.6 


1.30 ±0.03 


0.90 ±0.09 


2.20 ±0.10 


3C 147 


5.7 


1.26 ±0.03 


0.75 ±0.09 


2.01 ±0.09 


3C 147 


5.8 


1.26 ±0.04 


0.61 ±0.07 


1.87 ±0.08 


3C 147 


5.9 


1.21 ±0.04 


0.76 ±0.08 


1.97 ±0.09 


3C 147 


6.0 


1.20 ±0.03 


0.75 ±0.07 


1.95 ±0.08 


3C 147 


6.1 


1.16 ±0.03 


0.60 ±0.07 


1.76 ±0.07 


3C 147 


6.3 


1.16 ±0.03 


0.52 ±0.07 


1.68 ±0.08 


3C 147 


6.4 


1.18 ±0.04 


0.54 ±0.09 


1.72 ±0.09 


3C 147 


6.7 


1.13 ±0.04 


0.73 ±0.15 


1.87±0.15 


3C 147 


7.0 


1.04 ±0.04 


0.51 ±0.07 


1.55 ±0.09 


3C 147 
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Table 5 

Flux density measurements of Arp 220 



V (GHz) 


Flux (Jy) 


Calibrator 


1.1 


0.36 ±0.02 


3C 48 




0.34 ±0.01 


3C 286 


1.2 


0.37 ±0.02 


3C 48 




0.33 ±0.01 


3C 286 


1.4 


0.32 ±0.01 


3C 286 


1.5 


0.26 ±0.02 


3C 286 


1.6 


0.31 ±0.03 


3C 286 


1.7 


0.30 ±0.01 


3C 286 


1.8 


0.31 ±0.01 


3C 286 


1.9 


0.30 ±0.01 


3C 286 


2.0 


0.30 ±0.01 


3C 286 


2.5 


0.27 ±0.02 


3C 48 




0.29 ±0.01 


3C 286 


2.6 


0.28 ±0.02 


3C 48 




0.28 ±0.01 


3C 286 


2.9 


0.27 ±0.01 


3C 286 


3.2 


0.27±0.01 


3C 286 


3.3 


0.27±0.01 


3C 286 


3.4 


0.28 ±0.01 


3C 286 


3.5 


0.27±0.01 


3C 286 


3.6 


0.24 ±0.01 


3C 286 


4.3 


0.25 ± 0.02 


3C 286 


4.4 


0.23 ±0.02 


3C 286 


4.0 


n 0/1 -u n no 


0\u ZOO 


4.6 


0.24 ±0.02 


3C 286 


4.7 


0.26 ±0.02 


3C 286 


5.0 


0.24 ±0.02 


3C 48 




0.23 ±0.01 


3C 286 


5.2 


0.23 ±0.01 


3C 286 


5.4 


0.24 ±0.01 


3C 286 


5.7 


0.20 ±0.02 


3C 48 




0.22 ±0.01 


3C 286 


6.0 


0.22 ±0.01 


3C 286 


6.2 


0.29 ±0.03 


3C 286 


6.3 


0.20 ±0.03 


3C 48 




0.20 ±0.01 


3C 286 


6.7 


0.23 ±0.03 


3C 286 


7.0 


0.19 ±0.02 


3C 286 
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Table 6 

High-frequency observations. See for discussion. 



1/ Flux Telescope Tel. Kind^ dhAS^ Reference 

(GHz) (Jy) (arcsec) 



M82 



22.4 


1.00 ±0.15 


VLA D-cfg. 


I 


45 


Carlstrom & Kronberg (1991) 


23 


1.36 ±0.07 


Effelsberg lOOm 


SD 




Kronbere et al. ( 1979) 


24.5 


1.19 ±0.03 


Effelsberg lOOm 


SD 




Klein et al. ( 1988) 


32 


1.02 ±0.06 


Effelsberg lOOni 


SD 




Klein et al. (1988) 


42.2 


0.82 ±0.16 


VLA C-cfg. 


I 


43 


Rodn'euez-Rico et al. (2004) 


87.2 


0.51 ±0.08= 


Kitt Peak Um 


SD 




Jura et al. ( 1978) 


92 


0.59 ±0.09 


BIMA 


I 


50 


Carlstrom & Kronberg (1991) 


92 


0.67 ±0.10 


OVRO mm Arr. 


I 


55 


Seaauist et al. (1996) 



NGC 253 Core Component 



22 


0.49 ±0.05 


ATCA EW367-cfg. 


I 


61 


Ricci et al. (2006) 


23.7 


0.56 ±0.06 


VLA DnC-cfg. 


I 


60 


Takano et al. (2005) 


24 


0.52 ±0.05 


ATCA EW367-cfg. 


I 


56 


Ott et al. (2005) 


94 


0.32 ± 0.03 


BIMA A/B/C-cfg. 


I 


103 


Pentr et al. ( 1996) 



Arp 220 



22.5 


0.090 ±0.006 


VLA A-cfg. 


I 


2 


Zhao et al. (1996) 


23.1 


0.034 ± 0.007^^ 


VLA A-cfg. 


I 


2 


Norris (1988) 


42.2 


0.044 ±0.004 


VLA C-cfg. 


I 


43 


Rodrieuez-Rico et al. (2005) 


87 


0.030 ± 0.006^^ 


RAINBOW AB-cfg. 


I 


8 


Imanishi et al. (2007) 


87 


0.035 ± 0.007^^ 


IRAM PdBI 


I 


40 


Radford et al. (1991) 


97.2 


0.061 ±0.010 


IRAM 30m 


SD 




Anantharamaiah et al. (2000 



''The letter "I" denotes an interferometer, "SD" denotes a single dish. 
''Largest angular scale to which the observation is sensitive. 



'^Converted to flux density scale of lBaars et al\ (|1977l ). 
'^Uncertainty of 20% adopted. 
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Table 7 

Results of spectral modeling. 



Parameter^ 


M82 


NGC 253 


Arp 220 


\J "X.T 


1.21 


1.02 


1.17 


A 
B 
C 
n 


0.93 ±0.01 
-0.56 ±0.02 
-0.12 ±0.03 
04+° °^ 


0.50 ±0.01 
-0.51 ±0.03 
-0.16 ±0.03 

0.27 ±0.06 


n ^0+0.007 

n ni+0.03 
U.Ul_o.o4 

_r. O1+0.02 

"••'■'--0.04 

< 0.06 


as 
A 

/th 


-0.73 ±0.04 
-0.24 ±0.07 
0.06 ±0.02 


-0.74 ±0.04 
-0.33 ±0.06 
0.07t°:°? 


^•^^-0.05 

-0 62+0-04 

"•"^-0.08 

< 0.02 



"^Uncertainties in the parameters are marginalized la values determined via Monte Carlo simulations with 
10000 realizations for each source. 
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4 




100 200 300 400 500 600 700 800 900 1000 

Spectral Channel 

Fig. 1. — An assessment of the RFI present during a block of 1.7 GHz observations. This plot is typical 
of those produced during execution of the semi-automated RFI excision algorithm described in Solid 
lines show summed amplitudes of the 1.7 GHz data obtained during segment #14 as a function of correlator 
channel, corresponding to frequencies ranging from 1.65 to 1.75 GHz. Shaded regions indicate RFI-afflicted 
regions detected by the peak-finding algorithm. Plots such as these were reviewed and sometimes edited to 
improve the flagging {e.g., to widen the flagged region around the broad peak centered near channel 280). 
Channel 512 corresponds to a sky frequency of 1.700 GHz and the channel width is 102.4 kHz. The sinusoidal 
features across the bandpass are due to the ATA digital filter. 
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10^ F 




u 



I I I I 



y 



1 1 1 1 




1 1 1 



_I_L 



3 3.1 3.2 3.3 3.4 3.5 3.6 3.7 3.8 

log^o Frequency (MHz) 

Fig. 2. — The RFI environment of the ATA, in schematic form. Solid lines show a processed version of the 
summed data amphtudes, as described in ^ 33. II The raw 100 kHz channelized amphtudes were processed 
to increase clarity while accurately portraying narrow RFI peaks. Gaps in the data are due to a lack of 
observations, usually because the RFI was known to dominate astronomical signals at certain frequencies. 
The shaded regions indicate bands afflicted with especially strong and persistent RFI. 
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Fig. 3. — Average data retention rates over all observations, as a function of observing frequency. Entirely 
discarded datasets are not shown. Triangles denote the fraction of antpols used during imaging as compared 
to the number present in the raw data. Circles denote the same for individual visibilities. The shaded regions 
represent bands afflicted with strong, persistent RFI. The maximum possible visibility retention rate is ~75% 
due to the flagging of edge channels, shadowed antennas, and short baselines described in ^ 33. II 
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0''52" 50" 48" 46" 44" o''52" 50° 4B" 46" 44" 

HA (J2000) RA {J2000) 



Fig. 4. — ATA maps of NGC 253 at 1.3 GHz. The left panel shows the reference image for this source and 
frequency while the right panel shows an image generated from a 70 s integration made during Segment #3. 
The total integration time for the reference image is 250 s. Contour levels in both panels are 30 mJy x 2"/^ 
for n =- 0, 1, 2, . . . , 13. The RMS noise levels of the reference image and the snapshot are 9 and 11 mJy, 
respectively. Although the longer integration time and more complete u-v coverage of the reference image 
lead to a noticeable increase in its quality, the total fluxes of NGC 253 measured in the two images agree to 
2%. 
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_l I , I , I I I I \ I I I I . I , ^1 I . I 1 I I I I I 1 I I 1 I I L 

lO'^ia'" 6"^ lO'^O'" 54" 48"" 42"" 10^12" lO^O'^ 54"" 4B" 42" 

RA (J2000) HA (J2000) 



Fig. 5. — ATA maps of M82 at 1.0, 1.5, 2.0, and 4.5 GHz {top left, top right, bottom left, and bottom 
right panels, respectively) . The depicted maps are the reference images at each frequency reimaged onto a 
common coordinate grid. Circles: half-power points of the ATA beam. Filled ellipses: synthesized beamsizes. 
Contour levels in all panels are 50 mJy x 2"/^ for n = 0, 1, 2, . . . , 13. The RMS noise levels of the images 
are 17, 13, 10, and 11 mJy. The number of sources seen decreases as frequency increases due to both the 
shrinking primary beam of the telescope and the general downward slope of the radio spectra of celestial 
sources. 
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200 




Fractional Self- Calibrated Residual (%) 



Fig. 6. — Histogram of fractional flux residuals between calibrator observations and analytic calibrator 
models. The fluxes for the calibrator observations are derived from u-v modeling of calibrator data that 
have been self-calibrated to the calibrator reference images. 
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0.01 0.02 0.03 0.04 0.05 0.06 

Relative Uncertainty r] 

Fig. 7. — Contours of = k in the rj-aa space of uncertainty augmentation parameters. See ^3A\ for 
discussion. Solid line: the contour for values derived from aU repeated flux measurements. Dashed line: 
the contour for values derived from measurements of sources with Si, > 3 Jy. Dot-dashed line: the contour 
for values derived from measurements of sources with 5*^ < 1 Jy. 
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Fig. 8. — Spectrum of M82 with model. The two panels show different views of the same data. Cir- 
cl es: new measureme n ts pre sented in this paper. Other symbols i ndicate data from the literature: 



Carlstrom &: Kronberg (1991) and citations therein: 



( 19781 ) : lKronberg et (■197a ): lLeeuw fc RobsonI (|2009[ ): lRodrfguez-Rico aLl (|20o4 ): ISeaquist eraLl (|l996l ). 



Klein et al. (Il988h and citati o ns therein; Jura et al. 



These measurements are made with various telescopes in various configurations and, due to the properties of 
interferometers, are not always directly comparable. Squares: data obtained with single dishes. Triangles: 
data obtained with interferometers. Filled symbols: data used in the fits described in ii4.1.1l Open symbols: 
data not included in the fits. Shaded region: the Scr uncertainty region of the model spectrum as described 
in §4.21 Dashed line: the model nonthermal contribution. Dash-dotted line: the model thermal contribution. 
Dash- double- dotted line: a power-law fit to the FIR data. Solid line: the sum of the thermal, nonthermal, 
and FIR power law contributions. See ^4.2.11 for further discussion of the model results. 
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Fig. 9. — Spectrum of t he core componen t of NGC 253 with mod e l. The sy r nbols are as in Figure El 
External data are f ound inlBeck al\ (|l979l):lHildebrand et ad (Il977l): lOtt et al\ (|2005f) : |Peng et al\{imm 
Ricci et al\ (|2006h : iRieke et al\ (|l973h : Ixakano et al\ (|2005l ): IWeifi et al\ (|2008h . See j |4X2l for discussion. 
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Fig. 10. — Spectra of the extended component and total flux of NGC 253. The mod el used in F i gure P 
is overlaid to eas e comparison. Externa.1 dat a are found in iBeck et aLl (ll979^: fG cldzahl er fc Witzell (Il98lf) : 



Jura et al\ (|l978l ): IWright fc Otrupcekl (|l990l ): I Strickland et al\ (|2004l ): IWeifi al\ JgOOi). Triangles: fluxes 
of the extended component of NGC 253. Circles: sum of extended and core component fluxes. Other symbols 
are as in the previous plot. 
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Fig. 11. — Spectrum of Arp 220 wit h model. The symbols are a s in Figure [H See |34.2.3 | for further 
discus si on. Externa l data, a re found in Anantharamaiah et 'al\ (2000[): Carico et al\ ( 1992 ): Dunne fc Bales 



(pOoT); ' Bales et all (Il989l): llmanishi a/.l (120071): iNorrid (Il988l): iRadford et a l ('l99l V iRigopoulou eToI 



1996,):,Rodrfguez-Rico et all (|2005l) : IScoville et all ^199]h : Isopp fc Alexander! p991.) : .Woody et all (|l989l ): 



Zhao et all (|l996f) . 
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Fig. 12. — Flux density ratios of calibrator observations and spectral models. See §3.61 for definitions. 
Dashed lines: ratios of the VLA 1999.2 models. Solid lines: ratios of the spectral models adopted in this 
work. Squares and diamonds: directly observed calibrator ratios. The good agreement between the original 
and inverted ratios indicates that there is negligible flux loss in the secondary calibration process f ii3.2.2p . 
Triangles: calibrator ratios computed indirectly from the ratio of two observations of a science target. Xs: 
measurements of Rf^ 286 ina-de at the VLA in 2008 September. The observations of Rf^^j 2se were made in 
2009 June and exhibit discrepant behavior; see §3.61 for discussion. 
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B 

Fig. 13. — Observationally and theoretically derived values of the spectral curvature parameters B and C. 
Contours: la (inner) and 2a (outer) confidence regions for the parameters as deri ved from the modeling 
described in 21 Shaded region: the range of B and C values taken by the model of T+061 over a range of 



parameters appropriate for M82 as described in ^ The equivalent region for the NGC 253 parameters is 
virtually identical to the one shown. 
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Fig. 14. — Values of 77 and h^ff for which spectral curvature parameters B and C derived from the T+OGI 
model are consistent with observations of M82 to within 2a. Upper shaded region: the allowed parameters 
for p — 2. Lower shaded region: the allowed parameters for p = 2.3. Solid line: the value of 77 such that 
B = i ^min. D ashed line: the value of 77 such that the synchrotron cooling timescale Tgyn and the escape time 
Tesc of T+Od are equal, assuming v = 1.43 GHz, /loff — 100 pc, and the gal actic wind speed = 500 km s~^. 
This defines the approximate lower limit in r/ of the validity of the lT+06l numerical model. Dot-dashed line: 
the value of 77 such that B = B^q. 
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Fig. 15. — Values of ij and hos for which spectral curvature parameters B and C derived from the T+OGI 
model are consistent with observations of NGC 253 to within 2 a. The symbols are as in the previous figure. 
Shaded region: the allowed parameters for p — 2. The lT+n6l model with p = 2.3 is excluded to a 99.99% 
confidence limit. 
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